Effects of semiclassical spiral fluctuations on hole dynamics 
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We investigate the dynamics of a single hole coupled to the spiral fluctuations related to the 
magnetic ground states of the antiferromagnetic Ji — J2 — J3 Heisenberg model on a square lattice. 
Using exact diagonalization on finite size clusters and the self consistent Born approximation in the 
thermodynamic limit we find, as a general feature, a strong reduction of the quasiparticle weight 
along the spiral phases of the magnetic phase diagram. For an important region of the Brillouin 
Zone the hole spectral functions are completely incoherent, whereas at low energies the spectral 
weight is redistributed on several irregular peaks. We find a characteristic value of the spiral pitch, 
Q = (0.7, 0.7)7r, for which the available phase space for hole scattering is maximum. We argue that 
this behavior is due to the non trivial interference of the magnon assisted and the free hopping 
mechanism for hole motion, characteristic of a hole coupled to semiclassical spiral fluctuations. 



I. INTRODUCTION 

The interplay between charge and spin degrees of free- 
dom in a 2D doped Mott insulator represents an impor- 
tant problem in condensed matter physicSfi>S. The strong 
constraint on the double occupancy gives rise to the non 
trivial coupling of the charge motion with the magnetic 
background. Within the context of the t — J model, 
Shraiman and Siggiai^ showed that hole motion produces 
a long range dipolar distortion of the staggered magneti- 
zation resulting in a spiral order characterized by a pitch 
proportional to hole doping. Further studies focused on 
the spiral stability under low doping^ and its possible 
relation with the spin glass behavior found in the su- 
perconducting cuprates. Regarding the charge degrees 
of freedom, it was suggested that already for the one 
hole case such a distorsion of the magnetic background 
would lead to an orthogonality catastrophe,^ signalled by 
the vanishing of the quasiparticle weight Zk. Subsequent 
worksf^ however, showed that such a dipolar distorsion 
could be compatible with the spin polaron picture where 
Zk 7^ for the whole Brillouin zone (BZ). 

In frustrated Mott insulators the single hole dynamics 
can change drastically. That is the case for the 120° Neel 
order of the triangular antiferromagnet (AF)^ where the 
spin polaron picture is not valid for an ample region of 
the BZ; while for the kagome AF, whose magnetic ground 
state is disordered, the hole spectral functions seem to be 
completely incoherent^ for the whole BZ. Therefore, the 
coupling of the hole with the underlying spin excitations 
is crucial. In the former case, the hole couples with the 
magnons above the 120° spiralling magnetization; while 
in the latter case it couples with an exponentially large 
number of singlets inside a triplet gap.^ 

A hole can also couple with more exotic magnetic 
excitations like deconfined spin-i spinous. The decon- 
fined scenario, based on effective fields theories, was orig- 
inally proposed in the presence of short range spiral 
correlation^*^ and, later, in the quantum critical point 
between Neel and valence bond solid phases Given 
that these magnetic correlations are present in the AF 



Ji — J2 — J3 Heisenberg model it is believed that it is an 
appropriated microscopic model for the search of decon- 
fined spinous.— Recently, it was hypothesized that by 
coupling a non magnetic vacancy to the magnetic fluc- 
tuations of the above disordered regimes, in particular 
the plaquette ones, the strong reduction of the vacancy 
quasiparticle weight would give a measure of the degree 
of spinon deconfincmcnt. = Even if the idea seems to be 
appealing the QP weight reduction can be interpreted 
within a more conventional scenario 1^ 

Very recently, the phase diagram of the J1 — J2 — J3 
model was studied exhaustively with functional renor- 
malization group, coupled cluster method, and series 
expansion^**. The good agreement among these comple- 
mentary techniques allowed to establish the location of 
the quantum disordered (QD) regime in an important 
region of the parameter space (see right panel of Fig. 
[T]). Interestingly, in this QD regime substantial plaque- 
tte and short range incommensurate spiral fluctuations 
were found, analogously to early exact diagonalization 
studies on flnite systems. For this reason, a non trivial 
hole dynamics is expected in the QD regime of the above 
model. 

Motivated by this issue we investigate, as a first step, 
the spectral function of a hole coupled to the spiral fluc- 
tuations of the Ji~ J2 — J3 model. Using the t—J model, 
solved with exact diagonalization and within the self con- 
sistent Born approximation (SCBA),^^'^^ we found that 
i) for some regimes the sole inclusion of spiral fluctua- 
tions in the SCBA is enough to describe quite well the 
exact hole spectral functions on finite systems and ii) 
that a non trivial hole dynamics is observed under the 
effect of semiclassical spiral fluctuations in the thermo- 
dynamic limit. In particular, we found that in the strong 
coupling regime { J/t < 1) the spin polaron picture is not 
valid for an ample region of the BZ owing to the non 
trivial coupling of the hole with the spiral magnetic fluc- 
tuations which have been considered within the linear 
spin wave (LSW) approximation. Notably, in the weak 
coupling regime, { J/t ~ 10), where only few magnons are 
involved, it is observed a strong reduction of the quasi- 
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particle weight averaged on the entire BZ, Zav < 0.6 (see 
Fig. [SJ. In general, the QP coherence is lost due to 
the strong interference of the magnon assisted and the 
free hopping mechanisms for hole motion of the effective 
Hamiltonian, which increases the available phase space 
for hole scattering. By fine tunning frustration, we found 
that this effect is maximum when the magnetic pitch is 
around Q = (O.Ttt, O.Ttt). In many cases this spiral pitch 
corresponds to quite well ordered spirals, so we suggest 
that the QP vanishing is a direct consequence of the hole 
coupling with semiclassical spiral fluctuations. 

The paper is organized as follows: in Sec. IIA we 
briefly resume the magnetic phase diagram of the Ji — 
J2 — J3 model along with the LSW calculation. In Sec. 
IIB we present the t — J model treated within the self 
consistent Born approximation. In Sec. Ill we compare 
the spectra predicted by the SCBA with exact diagonal- 
ization on finite clusters and analyze the spectra in the 
thermodynamic limit. In Sec IV we close with the con- 
clusions. 



II. MODELS AND METHODS 

A. Ji — J2 — J3 model 

The Ji— J2— J3 AF Heisenberg model on a square lattice 
is defined as: 

Hj = Ji^Si • Sj+ J2^ Si • Sfc+ J3 ^ Si • S„i, (I) 

<ij> <ik> <im> 

where <ij>, <ik>, and <im> indicate sums to first, 
second, and third neighbors, respectively. The classi- 
cal phase diagram of this model is shown in the left 
panel of Fig. [T]and contains Neel (7r,7r), spirals {Q,Q) 
and {Q,tt), and collinear (tt, 0) or (0,7r) phasesJ^ The 
thick and thin lines indicate continuous and discontinu- 
ous transitions, respectively. Actually, the collinear (tt, 0) 
and (0,7r) phases are degenerated with an infinite num- 
ber of interpenetrating Neel states but it is known that 
quantum fluctuations select the collinear ones by the well 
known order by disorder phenomenoufi^ 

For spin S=^, a semiclassical spin wave study predicts 
an intermediate disordered phase along the classical tran- 
sition lines between the classical phases^i^ During a long 
time the quantum S = case of the Ji — J2 model has 
been intensively studied in the literature^ Only recently, 
however, the complete quantum phase diagram has been 
investigated with more sophisticated techniques. In the 
right panel of Fig. [T] it is sketched the quantum phase 
diagram of the model for 5 = i, recently foundji^ In 
particular, in the quantum disordered regime there is ev- 
idence of short range plaquette and incommensurate spi- 
ral correlations. As we have mentioned in the introduc- 
tion in the present work we will consider the effects of the 
spiral fluctuations on the single hole dynamics and leave 
the effects of plaquette fluctuations for a future work. 
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FIG. 1: Phase diagram of the J1—J2—J3 model. Left panel 
is the classical one with the magnetic phases characterized 
by the spiral pitch indicated in parenthesiaiS.. Right panel: 
sketch of the quantum phase diagram for S—^ recently found 
in Ref.ii 

Since the hole dynamics depends mostly on the struc- 
ture of the magnetic fluctuations, rather than on their 
long - or short- range character, and given that the spiral 
pitchs practically coincides with the classical ones,— we 
assume the magnetic background within the linear spin 
wave theory. For spiral phases the spin operators of Hj 
([T|) are expressed with respect to a local axis pointing in 
the classical direction of the spin at each site. Using the 
Holstein-Primakoff transformation in Eq. ([T]), it results 
the following quadratic Hamiltonian 

H.J = 7qaqaq + ^ ^1 ("q"-q + «-q«q) +^<^i (2) 
q q 

where 

7q — 2s Jg[cos^ cosq.J — cosQ.(5] 

<5>o 

= -2s Jg [sin^ cos q.S] , 

S>Q 

being Q the spiral pitch that minimizes Jk — Jge^^'^ 
with the sums on S extending up to third neighbors. 
After Bogoliubov transformation Oq — Uqdq + WqSLq, 
such as Uq— [ ^2^"" ]^' ^q = ~sign(/3q)[2^^^^]l, the 
Hamiltonian is diagonahzed as 

Hj ^Y Wqttljdq + ^ '^l + + ~)^'='- 

q q 

with a magnon relation dispersion 

ojI = s^[Jq - Jql[i(./q+q 4- Jq-^q) - Jq]. = 7q " Pq- (3) 

While the first member is the usual expression used in 
the literature for the spin wave dispersion,!^ t]^g second 



3 



one, in terms of functions 7q and /3q allows us to write 
down a more compact expression for the hole-magnon 
vertex interaction (see below). Notice that besides the 
three Goldstone modes at q = (0,0), and Q — zL{Q,Q), 
corresponding to the complete SO{3) symmetry rupture, 
at the linear spin wave level there appear two extra zero 
modes at k = (— Q, Q) and {Q, —Q) that reflect the lat- 
tice symmetry in the spectrum ([3]). In fact, a classi- 
cal spiral (Q, Q) is related to (Q, —Q) and (— Q, Q) by 
a global rotation combined with a reflexion about y and 
X axis, respectively.^^ Higher orders beyond linear spin 
wave theory will lift these degeneracies j22 nevertheless for 
our purpose it is enough to keep up to quadratic order 
since the main effects on the hole dynamics will be re- 
lated to the non-collinearity of the magnetic background 
signalled by the spiral pitch Q (see below). 



B. t-J model and SCBA 

The dynamics of a hole coupled with the magnetic ex- 
citations of a Mott insulator is properly described by the 
t — J model^ defined as 



of frustration leads to the vanishing of Zk. Regarding the 
hole-magnon vertex it is interesting to note that while 
for the Neel phase Mkq vanishes at q = (0, 0) and (tt, tt), 
for spiral phases it behaves as 



M, 



kq 



Af, 



k,Q+q 



-2j 



2i 



70 



R>0 



tji (q • R) sin 



Q.R 



7Q 



2w, 



Q+q 



cos k • R 



tR sin Q.R sin(k — ■^) • R 



i?>0 



for small q. As Wq~|q| and ajq+q'^jql, the above ver- 
tices behave as Afkq c>c |q| 2 and Mk,Q+q oc |q|~2 around 
(0, 0) and (Q, Q), respectively. To investigate the hole dy- 
namics we compute the hole spectral function Ak(w) = 
-(l/7r)ImG';^(w), where Gjj(w) =< AF|/ik[l/(w + i?7+ - 
Hef f)]h[.\AF > is the retarded hole Green function, and 
\AF > is the undoped magnetic ground state in the LSW 
approximation. In the SCBA, the hole self-energy is 
given by the following self-consistent equation^ii^iii 



Hf + H, 



<ij>(7 



h.c) + H.J (4) 



where the electronic operators are the projected ones, 
Cicr = (1 ~ ni-„)cia, that obey the no double occupancy 
constraint, and in our present case Hj is Eq. ([T]). To 
take care of the constraint we use the spinless fermion 
representationi^ for Ht , leading to the following effective 
Hamiltonian for the hole. 



Ek(t^) 



V- 

^ UJ - 



lAf, 



kql 



^k-q 



[uj - Wk-q) ' 



which must be solved numerically. The QP spectral 
weight can be calculated as Zk = [1 — 9Sk(w)/9w]~^|_Ek, 
where the QP energy is given by the equation i?k = 
ek + i?eSk(£;k). 



III. RESULTS: SCBA AND LANCZOS 



HeSf = ^fe^i^'' + 7^ 5](Afk,q/iLq/^k< + h.c.){b) 

k k,q 

The first term of Eq. ([5]) describes the free hop- 
ping of the hole without disturbing the magnetic 
background, and is characterized by the hole dis- 
persion ek = 2s ^jj^g 2tHC0s(^^) cos(k • i?) with R = 
(1, 0), (0, 1), in units of lattice space a. The second term 
describes the magnon assisted mechanism for the hole 
motion, and is characterized by the hole-magnon vertex 

A/kq = « (?7k-qWq " '7kWq) , 

where r]\^ = 2s '}2r >o 2iii sin(^p) sin(k • R). 

For a Neel phase the free hopping in Eq. ([5]) vanishes 
and it is widely accepted the spin polaron picture,— 1^ 
where the QP excitations are identified with the coherent 
propagation of the hole with the magnetic disturbance 
(magnon cloud). For spiral phases, instead, the free hop- 
ping processes interfere with the magnon assisted ones. 
In particular, we will show that the interference depends 
strongly on the spiral pitch Q, which for certain values 



A. Finite systems 

In this subsection we compare the hole spectral 
functions predicted by the SCBA and exact diagonal- 
ization on a cluster size of iV = 20 sites. For coUinear 
phases, (7r,7r) and (7r,0), we have already confirmed in 
a previous work a very good agreement between both 
techniques for the hole dynamics in the context of the 
Ji — J2 model. Here we are interested in frustration 
values that induce spiral correlations. In previous exact 
studies of the Ji~ J2— J3 modeli^ it has been shown a 
coexistence of spiral and plaquette fiuctuations which, 
in principle, would couple with the hole. Therefore, a 
careful comparison of the exact hole spectral functions 
with that of the SCBA would allow us to discern the 
relevance, or not, of the spiral fluctuations. It is worth 
to stress that for finite systems the parameter space is 
restricted to frustration values that give rise to spiral 
pitchs coincident with the momenta of the TV = 20 BZ. 
Even in this case, the strong finite size quantum fluc- 
tuations renormalize the classical magnetic wave vector 
Qc;. So, in order to perform a faithful comparison, we 
have assumed the position q™ of the maximum exact 
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FIG. 2: Lanczos (continuous) vs. SCBA (dashed) hole spectral functions for a 20 sites cluster, momentum khoie = (0, vr) (left 
panels) and khoie = (0,0) (middle panels), and J2/J1 = 0.2, Ji/t = 0.4. Upper (lower) panels correspond to J3/J1 — 1.5 
{J3/J1 = 0.45). Right panels: Exact static structure factor S'(q). In the SCBA calculation, the magnetic wave vector Q 
corresponding to the highest value of S{q) in the exact calculation has been chosen. M — —(0.8, 0.4)7r, and A = (0.8, 0.6)7r. 



structure factor S{q) as the actual spiral pitch charac- 
terizing the magnetic background within the SCBA. In 
Fig. 2 it is shown the comparison of the hole spectral 
functions with Lanczos and SCBA in the strong coupling 
regime, Ji/t = 0.4, along with the exact structure factor 
for different values of J3 / Ji at constant values of J2 and 
hole momentum lihoie- It can be seen that the agreement 
is reasonably good. In particular, when the magnetic 
structure factor is characterized by only one main peak 
like that at A = (|7r, |7r) (see upper panels of Fig. 2) 
the agreement is quite good, in contrast to the case of 



two peaks (see A and X = (tt, tt) at lower panels of 
Fig. 2). The latter is related to the strong competence 
of spiral and Neel fluctuations for this regime; however 
in the SCBA only the main spiral peak at momentum 
A is considered. Given that the exact results incorpore 
all kind of magnetic fluctuations, our results suggest 
that to describe the main characteristics of the hole 
spectral function, at least in the regimes considered, it 
is enough to take into account the coupling of the hole 
with magnonic excitations above the spiral correlations. 



In addition, in the spiral regime there is a strong reduc- 
tion of the quasiparticle weight in the low energy sector 
of the above spectra. This issue will be studied in the next 
subsection for the thermodynamic limit. 

B. Thermodynamic limit 

We have computed the hole spectral functions for clus- 
ter sizes up to 1600 sites. In complete agreement with 
previous works we found that for the frustrated Neel 
(tt, tt) regime the hole spectral functions show a coher- 
ent low energy peak for the whole BZ, which is identified 
with a well defined quasipaticle excitation, i.e. spin po- 
laron. This is shown in the dashed line of figure |3] where 



we have plotted the A(k, co) for a hole with momentum 
k/ioie = (O.Stt, O.Stt) at point A of the magnetic phase di- 
agram (see of Fig. [5]). For this case the QP weight is 
around Zk ~ 0.2 which means a considerable contribu- 
tion of multimagnon processes to the spin polaron wave 
function. As frustration increases the magnon dispersion 
bandwidth decreases, allowing to the hole to emit and 
absorb magnons more easily.^'^ Then, there is an increas- 
ing contribution of the multimagnon processes in the spin 
polaron wave function accompanied by a reduction of the 
quasiparticle weight. Therefore, as quantum fiuctuation 
are enhanced by frustration the QP weight and the mag- 
netization decrease monotonically. It can be seen that 
the QP weight remains finite even when long range order 
is destroyed, what confirms the idea that the hole dy- 
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FIG. 3: Hole spectral functions obtained by SCBA for 1600 
sites, khoie ~ (0.8,0.8)7r, J2 — O, and Ji/t = 0.4. Dashed (red) 
line: J3/J1 =0.1 (Neel phase). Continuous line: J3/J1 =0.425 
(spiral phase, Q — (0.7, 0.7)7r). Inset: low energy sector is 
zoomed. Label A and B, correspond to the points indicated 
in the magnetic phase diagram of Fig. [5] 



namic depends mostly on the structure of the magnetic 
fluctuations rather than on their long, or short, range 
character. 

In the incommensurate spiral phase we have found that 
the hole spectral functions are completely incoherent for 
an important region of the BZ. This is shown in the solid 
line of Fig. Owhere it is plotted ^(k, w) at point B of the 
magnetic phase diagram (see Fig. [S]) for the same hhoie- 
In the inset of Fig. [3] it is shown the peculiar structure of 
the low energy sector which is composed by several irreg- 
ular peaks where no signal of QP excitation is present. 
In the SCBA this feature of the hole spectral function is 
typical for an incommensurate spiral correlations and is 
related to the fact that there are two mechanisms for hole 
motion -the magnon assisted and the free one- whose in- 
terference may increase the available phase space for hole 
scattering, leading to the lost of QP excitations. In fact, 
this effect depends strongly on the spiral pitch which can 
be fine tuned by frustration. For instance, we have found 
that when the spiral pitch is around Q = (0.7, 0.7)n 
(case B of Fig. [3]) the effect is maximum with a van- 
ishing of the QP weight in approximately 50% of the BZ. 
It is worth to stress that here the local magnetization 
of the spiral state is about m ~ 0.24 so the non triv- 
ial hole dynamics is already present at a semiclassical 
level. We have found a similar behavior for a hole in- 
jected in other non-collinear magnetic backgrounds such 
as the 120° Neel order^ and canted Neel phases^i al- 
though the incommensurate spiral correlations seems to 
be more effective in destroying the QP excitations. It is 
worth noting that the strong reduction of the QP weight 
along with the rapid redistribution of the spectral weight 
on several mulipoles has been observed previously on fi- 
nite systems. Alternatively, these features were related 
to a scenario of spinon deconfinementiS which, of course, 




FIG. 4: QP weight averaged on the entire BZ, Zav = 
(1/A'^) jZk, as a function of J3/J1, predicted by the SCBA 
for 1600 sites in the weak coupling regime, Ji/t = 10. In be- 
tween the vertical lines m = within the LSW approximation. 
The dots B and C indicate the values of frustration for which 
Q « (0.7, 0.7)7r (they are displayed in the magnetic phase di- 
agram of Fig. O Inset: QP weight for \ihoie = (0.3, 0.3)7r 
computed with the SCBA (solid line) and Eq. ((6]) (dashed 
line) . 

is out of the scope of our present approximation. 

We have further investigated the effect of spiral fluc- 
tuations for the very weak coupling regime {Ji/t ^ 1) 
For a Neel phase the only mechanism for hole motion 
is magnon assisted and, owing to its high energy cost, 
it is expected that the hole reaches a quasi-static regime 
where the QP weight approaches unity. For spiral phases, 
however, the hole can move via the free hopping term of 
Eq.(l5]), reducing the QP weight even with a very low av- 
erage number of magnons promoted. This can be seen 
from the approximated expression of Zk for the weak cou- 
pling regimeyi^i^ 



In the inset of Fig. U it is shown the good agreement 
between the above expression (dashed line) and the QP 
weight computed with the SCBA (solid line) for different 
values of frustration. In particular, for the selected hole 
momentum, is strongly reduced once the magnetic 
background reaches the spiral regime. In principle, one 
can attribute the reduced value of Zk for spiral phases 
to the vertex behavior near {Q,Q), Afk,Q+q |q| ^. 
Nonetheless, we have checked that the most important 
contribution comes from the small values assumed by the 
denominators ek— Ck-q— Wq. In order to quantify globally 
this behavior we have plotted in Fig. 2] the QP spectral 
weight averaged on the entire BZ, Za«=(l/iV) X^k -^k, for 
Ji/t = 10. Again, it can be observed that as J3/J1 in- 
creases Zav decreases notably when the magnetic back- 
ground reaches the spiral regime. In particular, for 
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FIG. 5: The shaded area indicates the region of the magnetic 
phase space where the averaged quasiparticle weight Zav<0.6 
in the weak coupUng regime, {Ji/t — 10). For the strong 
couphng regime, Ji/t — 0.4, the SCBA predicts the vanishing 
of 2k in some area of the BZ. A = (0,0.1), B = (0,0.425), and 
C = (0.2,0.32). 

J2I Ji < 0.5 it is always observed a single minimum at 
Q ~ (0.7, 0.7)7r, analogously to that found in the strong 
couphng regime (Fig. For J2/J1 = (continuous 

line in Fig. |4]), Q occurs in a region where m k, 0.24, 
while for J2/J1 = 0.4 (dashed line) Q is inside the region 
where m = 0. In the latter we have modified the LSW 
approximation by introducing a gap in the magnetic dis- 
persion in order to describe the short range spiral order 
properlyi^ 

Finally, in Fig. [5] we plotted the approximated area 
of the magnetic phase diagram where the QP weight is 
strongly reduced, that is, zav < 0.6 for the very weak 
coupling regime. The fact that the shaded area is located 
at the spiral region of the magnetic phase diagram is a 
clear indication of the crucial effect of the spiral fluctu- 
ation on the hole dynamics. Since this is independent 



of the range of the magnetic fluctuations we expect the 
same effect for a hole coupled to the short range spiral 
fluctuation of the actual quantum magnetic phase dia- 
gram (see right panel of Fig. [T|). 



IV. CONCLUSIONS 

We have investigated the effect of semiclassical spiral 
fluctuations on a single hole dynamics. Based on the 
magnetic phase diagram of the Jj-Jt-Js model and solving 
the hole Green function with exact diagonalization and 
within the SCBA, we found that for the weak {Ji/t ^ 1) 
and the strong {Ji/t < 1) coupling regime there is char- 
acteristic value of the spiral pitch -Q « (0.7, 0.7)7r- for 
which the available phase space for hole scattering is 
maximum. Notably, for the whole spiral regimes ((Q, Q) 
and (Q, tt)) of the model we found a strong reduction of 
the averaged QP weight zav- In particular, for some mo- 
menta the QP weight vanishes and the spectral weight at 
low energy is redistributed on several irregular peaks, or 
multipoles. Even if in our study the spirals have been de- 
scribed semiclassically we think that this effect should be 
also observed in the short range spirals phases recently 
found in the disordered regime of the quantum phase 
diagramJ^ Similar features have been found in finite size 
studies of the same model, although the strong reduc- 
tion of the QP weight have been attributed to the spinon 
deconfinement inherent of the plaquette fiuctuations.— 
Based on our results and given that the alternative sce- 
nario for spinon deconfinement is based on short range 
spiral correlations)^ it would be important to go beyond 
the semiclassical description and investigate, within the 
context of the present model, the hole dynamics on spi- 
rals treated in terms of spinon excitations.^^ Work in this 
direction is in progress. 

We thanks A. E. Feiguin for his valuable help and use- 
ful discussions. This work was supported by CONICET 
under Grant PIP2009 N 1948. 
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